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Abstract. The probability that a number in many naturally occurring tables 
of numerical data has first significant digit (i.e., first non-zero digit) d is pre- 
dicted by Benford's Law Prob (d) = log 10 ^1 + i J , d = 1, 2 . . . , 9. Illustrations 
of Benford's Law from both theoretical and real-life sources on both science and 
social science areas are shown in detail with some novel ideas and generalizations 
developed solely by the authors of this paper. Three tests, Chi-Square test, total 
variation distance, and maximum deviations are adopted to examine the fitness 
of the datasets to Benford's distribution. Finally, applications of Benford's Law 
are summarized and explored to reveal the power of this mathematical principle. 



1. Introduction 

The significant-digit law of statistical folklore is the empirical observation that in 
many naturally occurring tables of numerical data, the leading significant digits 
are not uniformly distributed as might be expected, but instead follow a particular 
logarithmic distribution. Back in 1881, the astronomer and mathematician Simon 
Newcomb published a 2-page article in the American Journal of Mathematics de- 
scribing his observation that books of logarithms in the library were dirtier in the 
beginning and progressively cleaner throughout [1]. He inferred that researchers, 
be them mathematicians, biologists, sociologists as well as physicists, were looking 
up numbers starting with 1 much more often than numbers beginning with 2, and 
numbers with first digit 2 more often than 3, and so on. This ingenious discovery led 
him to conclude that the probability that a number has first significant digit (i.e., 

first non-zero digit) d is Prob (d) = log 10 (l H — ) , d = 1, 2... 9. In particular, 



d / 

his conjecture stated that the first digit is 1 about 30.1% of the time, and is 9 only 
about 4.6% of the time. That digits are not equally likely comes as somewhat of a 
surprise, but to claim an exact law describing their distribution is indeed striking. 
Passed unnoticed, the proposed law was discovered again and supported by empir- 
ical evidence by the physicist Benford who analyzed the frequencies of significant 
digits from twenty different tabled including such diverse data as surface areas of 
335 rivers, specific heat of thousands of chemical compounds, and square-root tables 
[2]. And this First Digit Law is known as Benford's Law today. But in recognition 
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of Newcomb's discovery, we can call it Newcomb-Benford's Law. This law applies 
to stock prices [3] , number of hours billed to clients [4] or income tax [5] as well as 
mathematical series [6]. And the tremendous practical values of Benford's Law were 
neglected until recently many mathematicians began to focus on the applications of 
this amazing phenomenon such as the design of computers and analysis of roundoff 
errors [7, 14, 15], as well as a goodness-of-fit against Benford to detect fraud [8]. 

In this paper, we will be dealing with heuristic argument and distributional property 
of the Significant-digit law in section 2; checking how data from various sources fit 
Benford's Law in section 3; discussing the application of Benford's law in section 4. 



2. Mathematical Formulation 
2.1 Heuristic Argument 

Pietronero and his colleagues gave a general explanation for the origin of the Ben- 
ford's law in terms of multiplicative processes in 2001 [3]. Here, the explanation 
was amended slightly in such a way that it can be used to explain the Benford's 
distribution not specifically to base 10. It stated that many systems such as the 
stock market prices which is discussed later do not follow the dynamical description 
by a Brownian process: 

N(t + 1)=£ + N(t), 
but rather a multiplicative process: 

N(t + 1) = £JV(t) 

where £ is a stochastic variable. By a simple transformation, we get 

lniV(t+l) = lnf + lniV(t). 

If we consider In £ as the new stochastic variable, we recover a Brownian dynamics in 
a logarithmic space; here we mean that a random multiplicative process corresponds 
to a random additive process in logarithmic space. This implies that as t — ► oo, the 
distribution Prob (In N) approaches a uniform distribution. By transforming back 
to the linear space we have 

J Prob(lniV)d(lniV) = J Cd(\nN) = C J ^dN, 

where C is a constant. 

It should be noted that Prob (N) = — is not a proper probability distribution, as it 
diverges or put in another way J °° — dN is undefined. However, the physical laws 
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and human conventions usually impose maximums and minimums. The probability 
that the first significant digit of N is n in base b is given by the following expression: 



Prob (n) 



r + ^iV _ln^ 
It TfdN ^ b 



log 6 



1 

1 + - 

n 



for any integer n which is less than b. We can review Prob (n) = log h ( 1 H — I as a 

V n. 

generalized expression of Benford's law to arbitrary base b. 

2.2 The Significant-Digit Law and Some Consequences 
2.2.1 Significant-Digit Law 

The Significant-Digit Law is 



Prob(L»i =d u ...,D k 



log 



10 



k-i 



vi=l 



where Di, D 2 , ■ ■ ■ Dk are the first, second . . . /c'th digits respectively. 

For example, Prob (Di = 1,D 2 = 2, D 3 = 9) = log 10 (l + (129)- 1 ) ^ 0.00335, which 
means there is a probability of 0.00335 that the first three significant digits are 129 
in a sample of Benford's distribution [11]. 

Hill has proved in his papers: "Scale- Invariance implies Base-Invariance" , not vice 
versa [9]; "Base-Invariance implies Benford's Law [10]"; "The logarithmic distribu- 
tion is the unique continuous base-invariant distribution [9]". He has explained the 
Central-limit-like Theorem for Significant Digit by saying: "Roughly speaking, this 
law says that if probability distributions are selected at random and random samples 
are then taken from each of these distributions in any way so that the overall pro- 
cess is scale (or base) neutral, then the significant-digit frequencies of the combined 
sample will converge to the logarithmic distribution [11]". 

2.2.2 Distributional Properties of D^s 

Based on the Significant-Digit Law, we analyze the statistics and distributional 
properties of a Benford's Distribution. 

2.2.2.1 Mean and Variance for Dk 

In this subsection, we compute the numerical values of the means and variances of 
Dk's using these expressions: 



E(D k )=Y, nPmb (D k = n) 
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Var (D k ) = Y, n2 Prob iPk = n) - E(D k ) 2 , 

n=l 

and tabulate them below. 

Table 1. Mean and variance of D k for k = 1 to 7 



k 


E(D k ) Var (D k ) 


1 


3.44023696712 (5) 


6.0565126313757 (6.67) 


2 


4.18738970693 (4.5) 


8.2537786232732 (8.25) 


3 


4.46776565097 (4.5) 


8.2500943647286 (8.25) 


4 


4.49677537552 (4.5) 


8.2500009523513 (8.25) 


5 


4.49967753636 (4.5) 


8.2500000095245 (8.25) 


6 


4.49996775363 (4.5) 


8.2500000000953 (8.25) 


7 


4.49999677536 (4.5) 


8.2500000000016 (8.25) 



This shows that by the Significant-digit law the mean of D k is approaching 4.5 which 
is the mean if the distribution were uniform and the variance of D k is approaching 
8.25 which is the variance if the distribution were uniform. 



2.2.2.2 Histogram of D k for k 
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Fig 1. Histograms of D k for k = 1 



Fig 2. Histograms of D k for k = 2,3,4, 5 



With the histogram, we illustrate that the distribution of the nth significant digit 
approaches the uniform distribution rapidly as n — > oo. 



2.2.2.3 Total variation distance of D k with respect to a uniform distri- 
bution 
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By calculating the total variation distance d(D±, U) of 
Dk from a uniform distribution on {0, 1, 2 ... 9} using 
these expressions: 

9 



Table 2. Total Varia- 
tion Distance of Dk from 
a uniform distribution on 
{0,1,2,. ..,9}. 



71= 
1 9 

d(D k ,U) = -J2 



Prob (L>i 



n=0 



->"5 



Prob(D fc = n)- — 



for fc 



for k ^= 1. 



As shown in the table, we observe that the is 
converging to the uniform distribution geometrically. 



k 


d(D fc ,C0 


1 


0.26872666 


2 


0.04702863 


3 


0.00488356 


4 


0.00048858 


5 


0.00004886 


6 


0.00000489 


7 


0.00000049 


k — > oo 


d(£> fc , £/)-»() 



2.2.2.4 Correlation coefficient of Dj and Dj 

Here, by calculating the correlation coefficient of D{ and Dj, 

_ Cov (D l ,D j ) 
Pd > d i ^Var(A)Var(Z) i ) 

for < i < j. we intend to investigate the dependence of one digit on another digit. 
We conclude that the dependence among significant digits decreases as the distance 
(j —i) increases. 

Table 3. Correlation coefficient of Di and Dj 





2 


3 


4 


5 


1 


0.0560563 


0.0059126 


0.0005916 


0.0000591 


2 




0.0020566 


0.0002059 


0.0000205 


3 






0.0000228 


0.0000022 


4 








0.0000002 



3. Illustrations 

The Chi-Square test formula to check the sample's goodness of fit to Benford's 
distribution is given as below, 

x *(8°) _ £ "° g(l + ^- Pr t (P ^" )>2 - S, 



n=l 



log (1 + £) 



where S denotes the sample size. The critical values for 8 degree of freedom at 5% 
and 1% level of significance are respectively 15.51 and 20.09. 

Another test, which we have used to analyze the sample, is the total variation 
distance (denoted as d\ here) of the sample from the Benford's distribution on 
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{1,2,... ,9}: 



d 



1 



71=1 



Prob (Di = n) - log(l + 



1 



n 



The maximum of deviations from a uniform distribution for each digit could also be 
considered as a test to check the goodness of fit: 



d 



max — max 

Kn<9 



Prob (Di = n) = log I 1 + - 

n 



3.1 Physical Constants 

Many literatures on Benford's law cite the table of physical constants as an example 
to illustrate Benford's law [9, 10, 11]. It was only Burke and his colleagues who 
actually attempted to check whether the physical constants would match Benford's 
law [12]. 

However, they only chose the con- 
stants from an introductory physics 
text. The sample size is too small 
to be significant. Here the 183 
constants from the 1998 Committee 
On Data for Science and Technology 
recommended complete listing of 
the fundamental physical constants 
(http:/ /physics. nist.gov/constants ) 
were analyzed for the first time. 



Table 4. First Significant Digit Test for the physical constants. 



Digits 


Counts in 


Predicted by 


Frequency Observed 


The Difference 




Sample Frequency 


Benford's Law 


in Sample 




1 


63 


0.3010 


0.3443 


0.0432 


2 


37 


0.1761 


0.2022 


0.0261 


3 


18 


0.1249 


0.0984 


-0.0266 


4 


15 


0.0969 


0.0820 


-0.0149 


5 


15 


0.0792 


0.0820 


0.0028 


6 


13 


0.0669 


0.0710 


0.0041 


7 


7 


0.0580 


0.0383 


-0.0197 


8 


7 


0.0512 


0.0383 


-0.0129 


9 


8 


0.0458 


0.0437 


-0.0020 


Sample 


Size: 183 




Chi-Square Test: 5.206 


Total variation distance: 0.0762 


Maximum of deviations: 0.0432 



3.2 Stock Prices & One-day Returns on Stock Index 
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Fig 3. Histogram on fitness of physical constant 
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Interesting results on the frequency of first significant digit of one-day returns on 
stock index were found by Ley [13]. However, investigations on stock prices were not 
done. According to Pietronero, due to the multiplicative process, a stock's prices 
over a long period of time should conform to Benford's Law. Here, we look at all 
stock's prices on certain days rather than a single stock's prices on a series of days. 
We collected data on local stock market SGX main board from http:/ /www. sgx.com 



Table 5. Analysis of Stock Prices 





Sample Source 


Sample Size 


x 2 


di 




A 


One-trading-day 
(13/10/2003) 


548 


2.23 


0.0277 


0.01734 


B 


20-trading-Day 
(10/10/2003- 7/11/2003) 


11,015 


45.5 


0.0255 


0.01501 


C 


32-trading-Day 
(10/10/2003- 26/11/2003) 


17,214 


79.3 


0.0284 


0.01834 



The observed frequencies roughly agree with the theoretical frequency predicted by 
Benford's law, indicted by the relatively small values of the total variation distance 
(0.0277, 0.0255 and 0.0284) of the sample and the Benford's distribution. However, 
if we performed the usual chi-square test of goodness of fit, we would reject the null 
hypothesis for datasets B and C, as the chi-square values are much greater than 
15.51 for the critical chi-square test value of 8 degrees of freedom at 5% significant 
level or 20.09 at 1% significant level. It is because of the large number of observa- 
tions - that is, the classical acceptance region shrinks with the sample size, given a 
significant level. On the contrary, the Benford's hypothesis would not be rejected 
for dataset A. 



One interesting thing, which in 
three sets of data d max always oc- 
curs at n = 2, attracts our atten- 
tion. By plotting the histograms J? 
of three datasets against a Ben- g 
ford's distribution, it is found iJ 
that the distributions of stock 
price data are in quite agree- 
ment among themselves but a lit- 
tle different from Benford's distri- 
bution. 
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We attribute this to the fact that a simple multiplicative process cannot precisely 
reproduce the complicated stock process thus the distribution of stock prices data 
may follow a curve that is steeper than the log curve of Benford's Law. 

3.3 Some Mathematical Series 
3.3.1 Fibonacci Series 

Here we provide a simple explanation for the conformance of Fibonacci series to 
Benford's law. Fibonacci series is characterized by the recursion relation: a n +2 = 
a n+i + a>n- This formula does not provide the explicit information of a general 
term from the series that may suggest the origin of Fibonacci series' conformance to 
Benford's law. However, it could be easily derived that say for a series with a± = 1, 
a 2 = 2: 

n-1 

/ 1 + V5 \ / 5 + 3\/5 \ / 1 - V5 ^ 

a n 




-3V5 
10 



l-v^ "' 1 / 5-3y / 5 \ 



Moreover, the magnitude of the second term is ( 1 O vo ) ( ) which is al- 
ways less than 1 and approaches zero when n gets larger and larger. Thus to analyze 
the first significant digit of a ra , we only need to analyze: 



a' = 



l + v^V 1 / 5 + 3^/5"' 



2 / V 10 / 

instead. Now we have lna' n = (n-1) In (j^^j + In ^ 5+ 1 3 v ^ ^ . When the sample 

size is large, distribution on n could be considered as uniform, with In ^ 1+ 2 V ^ ^ and 

m (^nr^) being constants, {lna^} will be uniform. It is then followed that the 
{a' n } and therefore {a n } will be uniformly distributed in logarithmic space and will 
conform to Benford's law. 



3.3.2 The Prime-number Series 

The prime number series is rather uniform below 100000, with the probability of 
each possible first significant digit being between 12.5% and 10.4%. Moreover, using 
the upper and lower bounds of function pi from the prime number theorem, it can 
be shown that the prime number sequence approximates a uniform distribution. 

3.3.3 Sequence a 11 in base b, where a G R, n G N 

The origin of Benford's distribution in a n is often attributed to its scale invariance 
which means that any power law relation is scale invariant: 

Prob(AiV) = /(A)Prob (TV). 
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An even simpler argument can be done. For a n in base b, where a is a constant, the 
sufficient condition for which the first significant digit of a n is d is d x b k < a n < 
(d + 1) x b k (that is \og b d + k < nlog b a < \ogb(d + 1 + k) for an integer k. Thus 
for each k, the probability for the first significant digit of a n being d is 

(log h (d+l) + fc)-(log 6 d + fc) ( d + l \ 

(k + 1) - k gn \ d J ' 

3.3.4 Factorial 

For the first 160 factorial numerical values, the chi-square test does not reject the 
Benford's distribution; however the total variation distance is too big for it to be 
considered as fit. 

3.3.5 Sequence of Power, i.e. n 2 , n 3 , n 4 , n 5 , . . . 

For n k series, as the constant k increases, the distribution becomes closer to the 
Benford's law. All values of the three tests get smaller as k increases, as shown in 
the analysis of n 2 ,n 5 ,n 20 and n 50 for n = 1 to 30,000. 

3.3.6 Numbers in Pascal Triangle 

We suspect that the numbers in Pascal Triangle obeys Benford's law, because it 
relates to Fibonacci Series closely while in another perspective it could be considered 
as a mixture of sequences of power. However in this paper, no data on Pascal triangle 
were analyzed. 



Tabic 6. Analysis of some mathematical series 



Sample 


Sample Size 


x 2 


di 


^max 


Fibonacci Series (7 series) 


10, 317 


0.0125 


3.8 x 10~ 4 


1.7 x 10~ 4 


Prime-number below 1,000 


168 


45.0 


0.2271 


0.1522 


Prime-number below 100,000 


9761 


3247 


0.4905 


0.1761 


1.007™, Specifically for n = 1 to 30,000 


30,000 


0.410 


1.2 x 10~ 3 


5.9 x 10~ 4 


1.007™, Specifically for n = 1 to 65,028 


65,028 


0.0329 


2.5 x 10~ 4 


1,2 x 10~ 4 


Factorial of 1 to 100 


100 


6.95 


0.0651 


0.04885 


Factorial of 1 to 130 


130 


8.97 


0.0871 


0.03492 


Factorial of 1 to 160 


160 


10.10 


0.0834 


0.03635 


n 2 , Specifically for n = 1 to 30,000 


30,000 


3.16 x 10 3 


0.1409 


0.09900 


n 5 , Specifically for n = 1 to 30,000 


30,000 


2.76 x 10 2 


0.0394 


0.03640 


n 20 , Specifically for n = 1 to 30,000 


30,000 


20.8 


0.0112 


0.00681 


n 50 , Specifically for n = 1 to 30,000 


30,000 


3.7 


0.0048 


0.00256 



3.4 Demographic Statistics and other social science data 
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The demographic data that the author used in this section are obtained from Energy 
Information Administration United States Department of Energy. 



Table 7. Analysis of demographic data 





Sample Size 


9 

X 


di 




GDP total for each 


1606 


33.0 


0.1009 


0.03352 


country in 1992-2001 










World Consumption of 


480 


40.5 


0.1113 


0.05936 


Primary Energy by Selected 










Country Groups, 1992-2001 










Population for each 


212 


6.59 


0.0792 


0.04160 


country in year 2001 










Crude Oil Production, 


4196 


32.2 


0.0383 


0.02166 


Import and Export, 










Stock Build Statistics for 










Each Country in 1990-2001 











The chi-square test rejects the null hypothesis (the data's obeys the Benford's dis- 
tribution) for all except the population statistics; while the total variation distance 
says that the oil statistics fit the Benford's distribution better than the other three. 

Other two sets of statistics in social science obtained from The World Affair Com- 
panion by Gerald Segal (1996) are also analyzed. The null hypothesis is accepted 
in both cases in the chi-square test. 



Table 8. Analysis of social science data 





Sample Size 


x 2 




^max 


Annual average number of 
people reported killed or 
affected by disasters per 
region and country in 1968-1992 


333 


11.05 


0.0715 


0.03222 


Weapon imports of 50 
leading recipients in 
US dollars in 1990-1994 


291 


11.06 


0.0832 


0.02872 



3.5 Numbers appeared in Magazines and Newspaper 

Among the data compiled in Benford's original paper [2], "numbers appeared in 
Reader's Digest" attracts our attention. Numbers appeared in a magazine or news- 
paper seems random at first glance; its conformance to Benford's law was deemed 
as a coincidence rather a rule. However, by referring to the Hill's paper, "selected 
at random and random samples . . . then the significant-digit frequencies of the com- 
bined sample will converge to the logarithmic distribution" [11], we hypothesize that 
numbers appeared in Reader's Digest could be considered as a mixture of data of 
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different distributions and more specialized magazines containing biased mixture of 
data will deviate from Benford's distribution more significantly. Here numbers ap- 
peared in both Innovation (science orientated) and The Economist (current affairs 
orientated, more on social side) were analyzed. It was found that neither conforms 
well to the Benford's distribution; nonetheless, the combination dataset conforms 
to Benford much better with total variation distance of 0.0474. This supports our 
hypothesis. 



Table 9. Analysis of numbers appeared in magazines 



Numbers appeared in 


Sample Size 


x 2 


di 


Innovation (Vol. 2 No. 4. 2002) 


152 


6.57 


0.0747 


The Economist (17th May 2003) 


449 


16.24 


0.0788 


Combination of the above two 


601 


9.22 


0.0474 



3.6 Other Datasets 

Other kinds of datasets, which are usually quoted as example of Benford's law 
include: survival distribution [14], the magnitude of the gradient of an image & 
Laplacian Pyramid Code [15], radioactive half-lives of unhindered alpha decays [11]. 
Datasets compiled by Benford in his original paper include: Rivers area, popula- 
tion, newspapers, specific heat, molecular weights, atomic weights, Reader's Digest, 
X— Ray volts, American league, addresses, death rate ... [2] They provide a wide 
range of examples to illustrate Benford's Law; however, the certainty of their con- 
formance to Benford's law varies largely. 

4. Application of Benford's Law 

The analysis of first few significant digit frequencies provides us a potential frame- 
work to examine the accuracy and authenticity of data values in numerical data 
set. 

4.1 The Design of Computers and Analysis of Roundoff Errors 

Peter Schatte has determined that based on an assumption of Benford input, the 
computer design that minimizes expected storage space (among all computers with 
binary-power base) is base 8 [16], and other researchers have started exploring the 
use of logarithmic computers to speed up calculations [17]. It will be crucial for the 
development of computing science in the information age. 

4.2 A Goodness-of-fit against Benford to Detect Fraud 

A. Rose and J. Rose developed a VBA code for detecting fraud under Excel working 
environment [8]. But the program can only do a first digit test but not a first 
significant digit test (for example, the 1 in 0.150 will not be counted by the program), 
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and the program will encounter error upon reading text (not numbers). We improved 
the program to eliminate these limitations of the original program. 

However, we argue that the use of a goodness-of-fit test against Benford to detect 
fraud is dubious. It has many limitations. Firstly, this method to detect fraud will 
not be accurate if the data in the account has very close built-in maximums and 
minimums. Secondly, a one-off embezzlement or very few entries of amendments 
made in an accounting record will not be detected using Benford's Law. Moreover, 
the public or specifically here the fraudsters learn very quickly; after a few cases 
of successful detections of fraud using Benford's law being reported, they may take 
special cares upon changing the entries in the account table in order to conform 
to Benford's law. Therefore, we conclude that Benford's law can only serve as an 
initial strategy to detect fraud and it will become less useful later. 

4.3. Others 

Mark Nigrini mentioned in his paper [4] that Benford's law could be used to test 
the accuracy of measuring equipment. Hill in his paper suggested Benford's law as 
a test of reasonableness of forecasts of a proposed model [9] . 

5. Conclusions 

Among the large range of samples we have analyzed, by using Chi-square test and 
calculating total variation distance and maximum of deviations, Fibonacci series, 
Sequence a™, population and number in magazines fit Benford's distribution well. 
Fitness to Benford's law of other data such as physical constants and stock prices is 
controversial, and it requires further investigations. Two common explanations for 
the origin of Benford's distribution are multiplicative process and scale invariance. 
One simple multiplicative process cannot reproduce most of naturally occurring 
phenomena thus by virtue of multiplicative process it cannot explain the origin of 
Benford's distribution for all datasets that are fit. The understanding of the origin of 
scale invariance has been one of the fundamental tasks of modern statistical science. 
How systems with many interacting degrees of freedom can spontaneously organize 
into scale invariant states will be the focus of further research. 

It is generally believed that the more chaotic and diversified the probability distri- 
butions are, the better the overall data set fits the Benford's Law. Thus, we put 
forward a bold conjecture that Benford's Law is an act of nature and can be used 
to indicate the randomness of our world. Therefore, it can be related to the sec- 
ond thermodynamics law of entropy. This will be an interesting area for further 
investigation. 
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